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Abstract 

We study spherically symmetric dynamical horizons (SSDH) in 
spherically symmetric Einstein/matter spacetimes. We first deter¬ 
mine sufficient and necessary conditions for an initial data set for 
the gravitational and matter fields to satisfy the dynamical horizon 
condition in the spacetime development. The constraint equations re¬ 
duce to a single second order linear “master” equation, which leads 
to a systematic construction of all SSDH initial data sets with certain 
boundedness conditions. Turning from construction to existence, we 
find necessary and sufficient conditions for a given spherically sym¬ 
metric spacetime to contain a SSDH. 
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1 Introduction 

Since their introduction four years ago by Ashtekar and his collaborators |3j 
|3], dynamical horizons have been found to be very useful for the study of 
the dynamical formation of black holes. Among other things, they have been 
used to study gravitational wave fluxes and their influence on the areas of 
black holes as they develop, and they have been used in building analogues 
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of black hole thermodynamics for astrophysical systems evolving into black 

holes 5|. 

How extensively do dynamical horizons occur in solutions of Einstein’s 
equations? While there have been studies, both numerical and analytical, of 
certain families of explicit examples pj, no attempt has yet been made to 
determine the parameter space of solutions which admit them. We do this 
here for the special case of spherically symmetric solutions. 

We obtain two types of results, on the construction of spacetimes with 
dynamical horizons (by constructing dynamical horizon initial data), and on 
locating dynamical horizons in given spacetimes. 

We first consider the problem of constructing initial data satisfying both 
the dynamical horizon condition, and the Einstein constraint equations. Propo¬ 
sition [21 parametrizes those initial data on / x S 2 which are spherically sym¬ 
metric, which solve the Einstein constraint equations, and which serve as a 
dynamical horizon for their spacetime development. We find that the col¬ 
lection T> of such initial data is parametrized by two free functions (£, r) on 
an interval and a pair of real constants (t/o,2/i)- That is, for every choice of 
these two functions, we can generate an explicit, initial data set of this sort; 
conversely, every data set in T> can be generated in this way. For simplicity 
we consider only bounded potentials, but there are many examples of dy¬ 
namical horizons with unbounded potentials im whose global properties 
can be studied using similar techniques. 

Note that since spherically symmetric solutions of the vacuum Einstein 
equations contain no dynamical horizons 1 , we consider here the Einstein 
constraint equations with generic matter, represented by an energy density 
field p and a radial momentum density field £. 

Second, we consider the problem of finding dynamical horizon(s) in a 
given (smooth, globally hyperbolic) spacetime. Lemma [ 7 | and Proposition |S] 
give conditions on the stress-energy tensor which lead naturally to a general 
existence result for dynamical horizons in the spacetime. 

Our arguments depend crucially on two quite interesting and novel iden¬ 
tities: a constraint (JHJ) on two components of the stress-energy tensor, and 
a reformulation of the spherically symmetric constraint equations as a lin¬ 
ear second order ordinary differential equation © for the volume function 
y = r 3 . We expect these expressions will play an important role in extending 

1 By Birkhoff’s theorem the vacuum region is locally Schwarzschild, which has vanishing 

outer expansion only on the horizons, which are not spacelike. 
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our results to non-spherical spacetimes. 


2 Preliminaries 

We work here with four dimensional spacetimes ( M,g,T ) satisfying the Ein¬ 
stein gravitational field equations G lw = 87tg T /ju , where M is a spacetime 
manifold, g is a Lorentz signature metric invariant under an 50(3) action 
and with Einstein curvature tensor G^ = R, w — \RgT^ u is a stress-energy 
tensor held which may or may not be a functional of a specific set of held 
variables, and G is the Newton gravitational constant. In our spacetime 
results (Section 01 we require that (M,g,T) satisfy an energy condition; in 
our initial data results (Section 0J), such a condition is not needed. 

Definition 1 An embedded hypersurface E 3 in a spacetime (M,g,T), is a 
dynamical horizon for the spacetime if the following conditions are met: 

(DH1) E 3 is spacelike; 

(DH2) E 3 is foliated by marginally outer-trapped 2-surfaces, i.e., surfaces such 
that the expansion of the outer future null vector field £ = £ + vanishes, 

e, = e + = 0 ; 

(DH3) the marginally outer trapped 2-surfaces are inner trapped, 9 n — 6_ < 0, 
where n = £_ is the inward future null vector field normal 2 to the 2- 
surfaces. 

A 2-surface with null normal expansions satisfying (DH2) and (DH3) is said 
to be marginally trapped. A hypersurface satisfying just (DH2) and (DH3) 
is called a marginally trapped tube [2j; however here we are concerned only 
with the dynamical horizon case. Under spherical symmetry we may assume 
the 2-surfaces are isometric 2-spheres of positive radius. 

A spacetime may generally contain a number of dynamical horizons. How¬ 
ever, as shown in [2], the foliation of a given dynamical horizon by marginally 
trapped surfaces is unique. It also follows from the results in j2] that if a 
given spacetime satishes the null energy condition (i.e., T(£,£) > 0 for every 
null vector £), then no spacetime region admits a local foliation by dynamical 
horizons. 

Normalised by g(n,£) = —2. 
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3 Constructing Initial Data 

In this section, we derive an algorithm for the construction of spherically 
symmetric initial data sets which satisfy both the constraint equations and 
the dynamical horizon condition. In doing so we determine the freely speci¬ 
fiable parts of such initial data, showing that there is a bijective correspon¬ 
dence between this free data on the one hand, and spherically symmetric 
dynamical horizon solutions of the constraints on the other hand. We then 
describe some example solutions, and we discuss some physically motivated 
restrictions one might impose on these solutions. 

3.1 Construction Algorithm and the Free Data 

An initial data set for the Einstein equations consists of (E 3 , 7 , K, p, J), where 
E 3 * is a three dimensional manifold (the initial slice), 7 is a Riemannian metric 
(the initial metric), K is a symmetric tensor (the initial extrinsic curvature), p 
is a non-negative function (87TG times the initial matter energy density), and 
J is a one form field ( 87 tg times the initial momentum density). Assuming 
this data satisfies the Einstein constraint equations, 

2 p = A, - K ah K ab + (tr 7 W) 2 , (1) 

J a = V 6 A a6 -V a (tr 7 X), (2) 

where i ? 7 is the scalar curvature of 7 , it follows (assuming appropriate matter 
evolution equations) [TU] [H] [fj that there is a unique spacetime (E 3 x/, g , T ) 
which i) satisfies the Einstein gravitational and matter field equations, ii) 
induces the initial data (E 3 , 7 , K, p, J ) on the embedded hypersurface 3 E 3 x 
{ 0 }, iii) is globally hyperbolic, and iv) contains (up to diffeomorphism) all 
other spacetimes which satisfy conditions i)-iii). This spacetime is called the 
maximal spacetime development of (E 3 , 7 , K, p, J). These results allow us 
to focus on the simpler initial data constraint system, rather than confront 
the challenges of the full evolution equations. Note, however, that if we are 
given a pair of solutions of the constraint equations, it is generally not easy 
to determine whether or not their spacetime developments are diffeomorphic. 

3 in the sense that 7 is the pullback of g to E 3 x {0}, K is the second fundamental form 

corresponding to E 3 x {0}, p = 87rG T(e_s_,e±) for e± the future-pointing timelike normal 

to E 3 x {0}, and J = 87rG T(e±, ■) 
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If we restrict our attention to initial data sets which are spherically sym¬ 
metric, with E 3 diffeomorphic to the three dimensional annulus I x S' 2 , then 
we may write 


7 = dR 2 + r 2 (R)(d9 2 + sin 2 6d(p 2 ), (3) 

K = \p{R)^ + {t{R) -%p{R))dR 2 , (4) 


and J = £(i?) dR , where R G (i?o, -Ri), $ G (0,7r) and </> G (0, 27t)). (Note that 
the coordinate i? measures radial geodesic distance, and r(R) is the spherical 
radius function.) Hence the full set of (spherically symmetric) initial data is 
parametrized by five real valued functions: r(R),p(R),T(R),p(R) and £,(R). 
Note that r = tr 7 A” is the mean curvature of the spacetime hypersurface E, 
whilst /i is the mean curvature in E of the 2-spheres S 2 of constant radius r. 
The constraint equations m for spherical data take the form 

= - +1 ^r 2 {\p 2 -pT + p), (5) 

(iT 

= (2r-3/i)—(6) 

We wish to study spherically symmetric solutions of the constraint equations 
which serve as dynamical horizons for their spacetime developments. The 
dynamical horizon condition 6 + = 0 takes the form 


2 r 


d 2 r 

dp 

r — 

dR 


2 dv 

0 i = 9 + = p+-— = 0, (7) 

r dR 

where p is the trace of the extrinsic curvature K over the S 2 tangent 2-planes. 
Substituting for p in © and © gives two equations for and eliminating 
this common term gives 

P + £ = ( 8 ) 

We thus obtain a somewhat peculiar, quite explicit, condition on the matter 
held part of the initial data which must be satisfied if the data set is to 
correspond to a dynamical horizon. 

Before stating the main result on the existence and parameterization of 
initial data, we note another curious property: substituting © and © in 
a leads to a linear “master” equation for y := r(i?) 3 : 


d 2 y 
dR 2 


+ 4-i^o 


(9) 
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Proposition 2 Suppose we are given functions r, £ € (^(M) and a solution 
V = y(R ) o/0 on an interval I = (i?i, R 2 ) on which y is everywhere positive. 
Defining r = y 1 ^ 3 , p by 0 and p by 0 gives a solution of the spherical con¬ 
straint equations 00. Conversely, 00 hold if the given data (■ r , p, r, p, £) 
satisfy 00 and 0 on an interval where r is everywhere positive. 

Furthermore, let y = y{R) be the solution of 0 with initial conditions 

y(R o) = Vo > 0, -^y(Ro) = j/i > 0, (10) 

and denote by I = (R 1 ,R 2 ), R 0 £ I, the maximal interval on which the 
condition r > 0 holds. Let (y, r, p, r, p, £) be the corresponding solution of 0- 
0. //i?i > — oo then we may normalise R so R\ = 0 and then lim^ 0 r(R) — 
0 and 7 is singular at R.\; likewise if R 2 < oo then r(R) = 0 and 7 

is singular at R 2 . 

The solution (y,r, p,r, p,£) on the interval (R\,R 2 ) specifies a dynamical 
horizon if additionally j^y > 0 on (Ri,R 2 ). 

Proof: Direct computation verifies the equivalence of the spherically sym¬ 
metric constraint equations 00 to the linear equation ©, augmented by 
formulas © for p, and © for p. This proves the first part of the proposition. 

The continuity assumptions r, £ £ C'°(M) ensure that the initial value 
problem for the linear ODE © together with the the initial conditions m 
has a unique global solution y € C 2 (M) which is positive in a neighborhood of 
R 0 ] hence lim^^ y(R) exists and the solution y(R) extends to R < R\. By 
maximality of / it follows that either Ri = —00 or Ri > —00 and y(Rf) = 0. 
Assuming now that R\ > — 00 , we may without loss of generality set R\ = 0. 
Uniqueness of solutions of linear ODEs shows that y\ = ^y(O) > 0, since we 
require that r is nowhere vanishing on I. Expanding in a Taylor series about 
R = 0 gives y(R) = r 3 (R) = y\R + 0(R 2 ), where the error term is controlled 
by ||t||oo + ||^| |oo- Thus 7 = dR 2 + (R 2 ^ 3 y^ 3 + 0(i? 4 / 3 )) {d0 2 + sin 2 9 df> 2 ) 
= ((3/pi) 2 r 4 + 0(r 7 ))dr 2 + r 2 (d9 2 + sin 2 0d<f> 2 ), so the metric is singular at 
the symmetry centre r = 0. Thus although the radial equation © is satisfied 
globally, the physical solution does not extend to points where r = 0. This 
holds for R 2 as well as R.\. 

Finally, since 


9 _ 


p 


2d 4 d 4 d 

r dR r dR 3 dR ^ ^ 


( 11 ) 
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when 9 + = 0, the dynamical horizon condition (DH3) is equivalent to -^r > 0 
and > 0. Observe that (DH3) serves only to restrict the solution to the 
sub-interval (Ri,R 3 ) C (i?i,i? 2 ) on which -j^y > 0; without it the above 
solution could be extended beyond i? 2 , to regions where 6 L > 0 . ■ 

REMARK: The dynamical horizon metric behaviour 7 ~ c 2 r 4 dr 2 + r 2 dfl 2 
is similar to that of a constant mean curvature hypersurface near a point 
singularity [3]. However, (JBJ) shows the spacetime curvature has components 
(ie. G(e±, £)) which are unbounded as r —> 0 and we show in Corollary H3 that 
r = 0 is a true spacetime curvature singularity. 


3.2 Examples 

As outlined above, the process of producing spherically symmetric initial 
data sets which satisfy the constraints and the dynamical horizon conditions 
from a given pair of freely chosen functions is very straightforward. We now 
illustrate this process with some examples. 


Maximal Data with Co-Moving Matter: Setting £ = 0 corresponds to 
choosing the matter to be co-moving with the surface orthogonal observers, 
while setting r = 0 results in the mean curvature of the initial data set 
vanishing, so the data is maximal. With these choices, equation © becomes 
simply y" = 0, which has the general solution y(R) = yi(R — Ro ) + Vo- I 11 
this case, we may without loss of generality choose R 0 = 0 ,y 0 = 0, and 
y 1 = a 3 (for a > 0) so y(R) = a 3 R. It then follows that r(R) = aRs, and 
consequently the metric takes the form 

7 = dR 2 + a 2 i?idR 2 , (12) 


where dVt 2 indicates the round sphere metric. Calculating y from © and 
substituting into the formula for K , we have 


K 


2 

3 a 3 R 


dR 2 - 


—^dti 2 . 

3aRs 


(13) 


Calculating p from ©, we have p(R) = 1 2 . We note that for this example, 

the interval on which the solution is regular is I = (0, + 00 ). The solution is 
singular at R — 0, where r = 0. We also note that since ^ = a 3 > 0, the 
entire solution satisfies condition (DH3). 
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CMC Data with Co-Moving Matter: We again set £ = 0, but set r 

equal to a non zero constant. The general solution 4 to the equation for y is 
y(R) — — (1 — e~ Rr ) + yo, which results in 

1 = dR 2 + (^(l-e- TR )+y 0 )ldn 2 . (14) 

T 

The expression for K is a bit messy, but straightforward to obtain. The 
behavior of the function y(R) for these examples depends on the sign of r. 
Recalling that t / 0 and y\ are both positive, we find that if r is positive, then 
y(R) has a zero for some negative value of R, and for large R it approaches a 
constant. If r is negative, then for negative R the function y either hits zero 
for some finite value R,\ or it approaches zero as R —> — oo, while for large 
positive values of R the function y(R) is convex increasing and unbounded. 
Thus in these cases, either / = (—oo, +oo) or / = (R 1 , +oo). 

Since R is a radial coordinate, which measures distance from the sym¬ 
metry axis, to obtain a spherically symmetric dynamical horizon from this 
example, we need to restrict the interval / to (0,+oo). When we do this, 
we find that while the metric is not degenerate as R —> 0, the curvature is 
unbounded as one approaches the axis. 

Note that since || = yie~ Rr > 0 for any constant r, the (DH3) condition 
is satisfied for the full range of I. 

Flat 3-metric: If we set r(R) = R, then (£, 7 ) ~ (M 3 ,<5). To see which 
expressions for K and p and £ are compatible with a flat metric 7 = 6, we 
substitute y(R) = R 3, into Q, giving 

6 R + 3 tR 2 - |£R 3 = 0. (15) 

Assuming co-moving matter £ = 0 leads to r = —j- { . From CD we find 
that y — — so K — —RdQ 2 and the matter density is p = jp. The 
data is clearly singular at R = 0, and nowhere else. Without the co-moving 
assumption, we find that equation © implies that 

3R(2 + tR — \£R 2 ) = 0. (16) 

Choosing, say, £ = gives p = r = — ^ and /j = —2/R as before. 
One then readily constructs K , noting its singular behavior at the symmetry 
axis. 

4 In this case, since the solution is not linear in R, setting yo = 0 is a restriction. So we 
let yo and y\ be any pair of positive constants. 



For all cases of this example, since ^ = 1 > 0, condition (DH3) for a 
SSDH is satisfied everywhere. 

3.3 Physical Restrictions 

In the examples just discussed, we have made restrictions on the choice of 
the free data £ and r based on mathematical convenience. In this section, we 
consider physically motivated restrictions, and some of their consequences. 

Our first physical condition on a dynamical horizon £ is that there is at 
least one point on the symmetry axis, which we may normalise to R — 0, 

r(R = 0) = 0. (17) 

Second, we recall that the definition of dynamical horizon requires that the 
foliation 2-surfaces S'f { should be contracting in the inward null direction 
n = £_, but that this condition plays no direct role in the solution of the 
dynamical horizon equations Thus we consider explicitly the inner 

trapped condition 

< 0. (18) 

As shown by PP, the inner trapped condition in spherical symmetry is 
non-evolutionary , meaning that if it is satisfied on a Cauchy surface then it 
holds throughout the globally hyperbolic development. We can also show 
that m follows from the outgoing momentum condition 

f > 0. (19) 

Finally, we recall that the conceptual picture of a dynamical horizon has 
it inside an event horizon, with the area function r having a finite bound 
r < 2m(oo) where m(oo) is the final Bondi mass. The examples above show 
that this condition is not satisfied by all solutions of the dynamical horizon 
constraint equations ©,©• This motivates the black hole condition, that 
there is a constant r* > 0 with 


r* := sup r, (20) 

(o ,R*) 

where R* is maximal , in the following sense: 

Definition 3 Given the data (r, £, yi, R*) E C°(I) x C°(I) x (0, oo) x (0, oo] 
specified on the interval I = (0, R*), let the corresponding solution ( y , r, /i, r, p, £) 
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E(r, ?/i, i?*) be constructed as prescribed in Proposition with R 0 = 0, 
y 0 = 0. If this solution satisfies y > 0 and 6_ < 0 on I, then (r, £, 2 / 1 , R*) 
is called a DH data set. Such a DH data set is DH-extenclible if there is a 
DH data set (f,£,yi,R) such that R > R* and f|( 0R ) = r, £|(o,i?) = £ and 
yi — 2/i- A DH data set (r, £, 2 / 1 , R*) is maximal if it is not DH-extendible. 
If the inner trapped condition is not assumed then we have a generalised DH 
data set, and a DH data set is weakly DH-extendible if it is extendible in 
the class of generalised DH data. 

Proposition 4 Suppose £(r, £, 2 / 1 , R*) = ( y , r, / 1 , r, p, £) is a dynamical hori¬ 
zon constructed by Proposition [H on I = (0, R*), so in particular, S satisfies 
the axis condition m- 

(i) Suppose (r, £, 2 / 1 , R*) is a generalised DH data set, so the inner trapped 
condition m is not assumed a priori. If the outgoing momentum 
condition m is satisfied on a sub-interval (0, R ) then the dynamical 
horizon is inner trapped on the same sub-interval, and liiri™^ 6L < 0. 

(ii) Suppose the black hole JH) and inner trapped m conditions hold on 
(0, R*), where R* is maximal. Then r(R) is increasing on (0 ,R*), 
limRfR* -j^r = 0, and either R* < 00 and 9~(R*) = 0, or R* = 00 . 

Proof: (i) It follows from 9 + = 0 and n = I — 2e# that = —4^ logr = 
— lih logy- Proposition [21 shows that if 2/(0) = 0 then for small positive R 
we have 6L ~ —| y'(Q)/R < 0. The master equation @ and 9 + = 0 show 
that on E we have 

±9_ = \9 2 _-T9„-2i < f 9 2 _-t9-, (21) 

since £ > 0 by assumption. Introducing h = —1/9- and integrating gives 

h(R) < e T{R) fh(e) + ^J e~ T{s) ds^j , 

where T(R ) = t(s) ds and e > 0 is chosen small enough that h(e) < 
1. Now r G C°(M) so h(R) is bounded above uniformly for finite R and 
thus 9 _ = —1/h is negative as required; even more, 9 _ is locally uniformly 
negative, since J^ +1 r(s) ds is bounded, for all R. 
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(ii) The inner trapped condition 9_ < 0 and 6L = —4^1ogr shows that 
■j^r > 0. The master equation ® and r, £ G C°(M) show -f^y = 3r 2 ^r 
is bounded locally uniformly on M, so lim^yR* -j^r < oo while r is bounded 
away from 0. Choosing any C° extension (f, £) of (r, £) gives an R > R* and 
a generalised DH data set (f, £, yi, R) extending E, by continuity of y at R*. 
Continuity of j^y implies that E is DH-extendible, contradicting maximality 
of R*, unless lim^f#* f^r = 0. Finally, if R* < oo then continuity of j^y 
shows that 6L(R*) = 0. ■ 

Remark: Proposition 0 provides a natural picture of a dynamical horizon 
inside an event horizon, starting at the central axis and extending radial- 
outward and geodesically complete, with future endpoint at i + . However, 
this picture relies on the assumptions we have made, in particular the inner 
trapped condition m and the continuity of the free fields r, Examples of 
marginally trapped tubes (cf. JIE]) show that the condition 9 + = 0 dehnes 
a hypersurface which may become null and timelike. In these examples, the 
dynamical horizon metric is inextendible as a spacelike hypersurface but not 
complete. The next results examine the solutions near the singular points 0, 
R*. 

Lemma 5 Suppose E = E(r, £, yi, R*) is a spherically symmetric dynamical 
horizon with R* < oo. Let (E R ,E_ i_) be a frame along E which is parallel 
transported by the spacetime connection. Then there is a finite boost (with 
parameter a) to the E-adapted frame (e#, e±) at R = R*, 

£r = \(a-a~ 1 )E L + \(a +a~ l )E Rl 

e_L = i(a + a~ 1 )E J _ + i(a — a~ 1 )E R . 

Furthermore, E has spacelike radial unit tangent vector e R at R = R* and if 
liniR^* 9_ < 0 then E is DH-extendible. 

REMARK: The spacetime parallel transport provides a reference frame, uni¬ 
formly equivalent to any other construction of spacetime frame, and therefore 
suitable for verifying whether or not e R goes null at R*. 

Proof: The extrinsic curvature K a b is defined by K(X,Y') = g(X,\7ye±), 
so X £R e R = K(e R , e R )e± and X CR e± = K(e R ,e R )e R . The parallel frame 
(E r , E_ l) satisfies V eR E R = X eR E± = 0. The corresponding null frames 
£± = e± ± e R and L± = E± ± E R are related by a boost with parameter 
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a = a(R ), ie. £± = a ±l L ±, and we find ^loga = K(eR,eR) — r «- fx = 
r + 2^ logr. Integrating gives 



so the boost a(R ) is bounded while r > 0 and / r < oo. Now 9_(R*) < 0 and 
9_ = —4^ logr show that j^r(R*) > 0, so the dynamical horizon solution 
extends as before. ■ 

A similar argument can be used to show that the Einstein tensor is un¬ 
bounded at the axis. 

Corollary 6 Suppose £ is a spherically symmetric dynamical horizon in a 
spacetime M, which satisfies the axis condition (n> and fy t(s) ds < oo. 
Then the central axis point R = 0, r = r(R = 0) = 0 of £ lies on the 
singular set of the spacetime M. 

Proof: Equation (JHJ) shows that the component G(e±, £) is unbounded, and we 
must show this holds in all frames near R = 0. It suffices to consider frames 
L, N which are spacetime parallel along £. As in Lemma 0 the tangent 
frames n, £ and parallel frame N, L are related by a boost with parameter 
a = a(R ); i.e., £ = aL and n = cT 1 ^, where a = r 2 0( 1) as R —> 0 by (jSSj) - 
It follows from r~ 2 = | (G(£,n) + G(£, £)) = |( G(L,N) + a 2 G(L,L)) that 
G(L,L) = 0(r~ 6 ) is the required unbounded coefficient. Note also that the 
boost n = 0(r~ 2 )N so £ is tangent to the past null cone at the central axis 
point; see (JU], Figure 3.). ■ 


4 Spacetime Picture 

Let (M,g,T) be a spherically symmetric spacetime which satisfies the Ein¬ 
stein equations G^ = 87tG T^ v with stress-energy T^ u , and consider the ques¬ 
tion of finding a dynamical horizon in M. If one were to exist in (M,g,T), 
then the initial data on this dynamical horizon would necessarily satisfy the 
conditions discussed in Sectional including (0). Since (JHJ) depends on the 
adapted frame (ej_,e^), it is not easy to check directly from the spacetime 
fields. One can however state an equivalent slice-independent condition: 
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Lemma 7 Let W C M be a spacetime region with stress-energy tensor T 
satisfying the strict null energy condition (SNEC) 

T(L , L) > 0 for all radial null vectors L in W. (23) 


Then the following are equivalent: 


(i) There is a future null frame (£,n) with g(£,n) = — 2 satisfying the 
condition 

\r 2 G(n, £) = 47rG r 2 T(n, £) < 1. (24) 


(ii) There is a spacetime frame (E±, Er) satisfying J2J). 

Remark: The Lorentz-invariant condition (121 arises frequently: for exam¬ 
ple, in If) the sign of the C-function (0 equation (2.3)) 


\r 2 G(£,£) 

1 -±r 2 G(£,n) 


(25) 


(or more accurately, the sign of C -1 ) determines the causal character (space- 
like/null/timclike) of the marginally trapped tube (MTT). Another example 
is the evolution equation jH] 

5£R-n 2 (l-lG ( <,n)r=) (26) 


for the areal function r 2 (u, v) in double null coordinates ds 2 = — Q 2 (u, v ) du dv+ 
r 2 (dd 2 sin 2 (6 ) ) df 2 ). 

Proof: To show that (i) =>■ (ii), we check the effect of a Lorentz boost 
on T(e_i,£). The boosted frame L = a£, N = a _1 n, E± = |(L + N) = 
\(a£ + a~ l, n ) satishes 

2 T(E ± , L ) = T(N, L ) + a 2 T(L, L) ; 


hence it follows from the SNEC PI and the Lorentz invariance of (ED that 
there is a unique a > 0 such that T(E±, L ) = |(47TG r 2 ) -1 , which is equivalent 
to the condition Pj). The converse (ii) => (i) follows from the above boost 
relation, together with the definitions of p and £ in terms of T. ■ 

It follows from this Lemma that, if we are given a spacetime which satishes 
the SNEC, then inequality PI is a necessary condition for it to contain a 
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dynamical horizon. We now look to find sufficient conditions for existence. 
Suppose that in our given spacetime (. M,g,T ) there is a non-empty region 
W C M in which the SNEC condition and (1221) both hold. From the lemma, 
it follows that through every point in W there is a spacelike vector e r such 
that the adapted frame (e±,eR) satisfies (0). The integral curves of the 
vector field e r give a foliation of W by spacelike hypersurfaces which satisfy 
(0). It follows from the results of Section 0 that if the spacetime contains a 
dynamical horizon, it must coincide with a leaf of this foliation. 

To see that condition 0), whilst necessary for a given leaf to be a dy¬ 
namical horizon, is not sufficient 5 , we use the constraint equations with 0 
and the expression 9 + — p + for the future outer expansion 9 + to derive 
a differential equation for 6 + . Specifically, writing © in terms of 6 + and 
eliminating ^ with ©, we obtain 

Te + + (!» + -r)e + = i- P -v ( 27 ) 

Clearly if we choose a spacelikc hypersurface with adapted frame 
relative to which p + £ = \ (which Lemma 2 guarantees that we always can 
do) then the right hand side of (1271) vanishes, and we see that (1271) admits 
the solution 9 + = 0. It then follows from ODE uniqueness that the sign of 
9 + is fixed on each leaf. Hence if p 6 W and 9+(p) = 0 then the leaf through 
p is a dynamical horizon. Finally, it follows from the Raychaudhuri equation 

D £ 9+ = -§0+»'8ttG T(£,£) 

and from the strict null energy condition that D £ 9 + < 0. Presuming the 
spacetime to be smooth, we find that there are four possibilities: 

(a) The region W has 9 + > 0 everywhere, so that it is entirely untrapped, 
and there is no dynamical horizon. 

(b) The region W has 9 + < 0 everywhere, so that it is entirely trapped, 
and there is no dynamical horizon. 

5 We have shown in Section 0 that while 0) holds, for any choice of a pair of free 
functions, we can construct initial data so that the hypersurface with that data is a 
dynamical horizon for the spacetime development of that data. It does not , however, 
follow that any initial data set satisfying m is a dynamical horizon: we show here that 
the additional condition 0+(p) = 0 for at least one point is necessary. 
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(c) The region W contains a single leaf with 9 + = 0 , but 6L > 0 everywhere 
on this leaf. So there is no dynamical horizon. 

(d) The region W contains a single leaf with 9 + = 0, and on some (possibly 
proper) subset of this leaf, < 0. This subset is a dynamical horizon. 

In summary, we have shown the following: 

Proposition 8 Let (M, g, T) be a smooth spacetime which satisfies the strict 
null energy condition W- 

1. If condition | is satisfied nowhere in M, then the spacetime contains 
no dynamical horizons. 

2. If there exists a non-empty connected region W C M in which condition 
P7P holds, and if the outward null expansion function 9 + has constant 
sign (with no zeroes) then W contains no dynamical horizon and is 
everywhere trapping or non-trapping. 

3. If there exists a non-empty connected region W C M in which condition 

holds, and if the expansion 9 + takes on both positive and negative 
values in W, then W contains a unique spacelike hypersurface with 
9 + = 0. The data on this hypersurface satisfy ( 0 ), and it is (entirely or 
partly) a dynamical horizon if and only if 6L <0 on all of it or part of 
it. In addition, this hypersurface divides W into a trapping region to 
the future and a non-trapping region to the past. 

REMARK: It is worth noting that, unlike the situation for dynamical horizons 
in non-spherical spacetimes, the spherically symmetric construction given 
here does not depend on any choice of (spherically symmetric) time slicing. 
This leads to a much stronger uniqueness statement (Proposition 0) than is 
available in the general case. 


5 Conclusion 

While the set of spacetimes which are spherically symmetric is very special, 
it is clear from physical considerations that, so long as matter is present, 
dynamical horizons can form in them. We determine here the necessary con¬ 
ditions for spherically symmetric spacetime to contain a dynamical horizon, 
and sufficient conditions as well. We also demonstrate their uniqueness, and 
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show how to systematically construct initial data sets which serve as dy¬ 
namical horizons for their spacetime developments (with two free functions 
parametrizing the collection of all such data sets). 

The situation regarding dynamical horizons in non-spherically symmetric 
spacetimes, is of course considerably more complicated. Still, a full under¬ 
standing of the spherically symmetric case provides a good first step toward 
uncovering the properties of dynamical horizons more generally. 
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